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Abstract 

The unipolar and bipolar macroscopic quantum models derived recently for instance in the 
area of charge transport are considered in spatial one-dimensional whole space in the present 
paper. These models consist of nonlinear fourth-order parabolic equation for unipolar case 
or coupled nonlinear fourth-order parabolic system for bipolar case. We show for the first 
time the self-similarity property of the macroscopic quantum models in large time. Namely, 
we show that there exists a unique global strong solution with strictly positive density to the 
initial value problem of the macroscopic quantum models which tends to a self-similar wave 
(which is not the exact solution of the models) in large time at an algebraic time-decay rate. 

1 Introduction 

The quantum hydrodynamic (QHD) model for semiconductors is derived and studied recently in 
the modelings and simulations of semiconductor devices, where the effects of quantum mechanics 
arise. The basic observation concerning the quantum hydrodynamics is that the energy density 
consists of one additional new quantum correction term of the order 0(e) introduced first by 
Wigner [27j in 1932, and that the stress tensor contains also an additional quantum correction 
part [T] related to the quantum Bohm potential [;3] 

with observable /3 > the density, m the mass, and e the Planck constant. The quantum 
potential Q is responsible for producing the quantum behavior. Such possible relation was also 
implied in the original idea initialized by Madelung [25] to derive quantum fluid-type equations 



2 H.-L. Li, G.-J. Zhang, M. Zhang, C.C. Hao 

in terms of Madelung's transformation applied to wave functions of the Schrodinger equation 
of the pure state. Recently, the moment method is employed to derive quantum hydrodynamic 
equations for semiconductor device at nano-size based on the Wigner-Boltzmann (or quantum 
Liouville) equation, refer to [13j for details. For more important progress on the derivation of 
macroscopic quantum models in terms of the entropy minimizer principle, one can refer to the 
recent interesting works [TJ [T3l [T5l [T6] and the references therein. 

Starting with the quantum hydro dynamical models and performing the relaxation limit 
asymptotical analysis, the macroscopic quantum (Drift-Diffusion) model is derived rigorously pi] 
for the model of the unipolar carrier, the methods employed therein can be generalized to gen- 
eral bipolar carriers. For positive charge density, these models are indeed nonlinear fourth-order 
parabolic equation for the unipolar case or the coupled nonlinear fourth-order parabolic system 
for the bipolar case. We would like to mention that according to the recent result in [6], the 
model (II. 4p can also be viewed as a relaxation limiting equation of the quantum fluid model 
which can be derived by the nonlinear Schrodinger-Langevin equation, for which the rigorous 
short time existence of weak solutions is proven recently in |2Uj . 

We are interested in the long time asymptotical behavior of solutions to the macroscopic 
quantum models in the present paper in the one-dimensional real line, and we shall show that 
the global classical solutions to the IVP (fOjl - ffO]) and the IVP pril) - p3]) admit the character 
of self-similarity in large time. In general, the typical bipolar macroscopic quantum model widely 
used in semiconductor modeling in one dimension is the following coupled nonlinear parabolic 
system 

dtPi - p{piU + (pi ( ^"^"^ ) ) +{-iy^HpiE), = 0, i = a,b, 

\ \ Vp^ JJ^ (1.2) 

^ = Pa - Pb, t>0, X £R, 
together with the initial data 

Pi{x, 0) = /Oi,o(x) > 0, X G M, pj,o(±oo) = p± > 0, i = a,b, (1.3) 

where pa, Pb > denote the macroscopic densities for electron and hole respectively [13], p{pi) 
is the pressure function depending on the density pj, e > is the scaled Planck constant, and 
E denotes the self-consistent electric field. We also use the symbols (—1)" = —1 and (—1)^ = 1 
for the simplicity of statements. 

In the absence of the electric filed (or in the so-called quasi-neutral domain), the initial 
value (IVP) problem ()1.2p reduces to the IVP problem for the following unipolar macroscopic 
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quantum model 



(1.4) 



p{x,0) = po{x) > 0, X eR, /5o(±oo) = p± > 0, 



(1.5) 



where p = p{x, i) > is the density of electron or hole and p = p{p) is the pressure function 
depending on p. It should be noted that a similar model (the DDLS model), which takes the 
form of (jl.4p but without the pressure term (i.e., p = 0), also arises in the study of interface 
fluctuations in spin systems, for instance [8j. 

There are recently many analysis results on macroscopic quantum models of the fourth-order 
parabolic type p.2|) or (|1.4|) and related models. For the Eq. (|1.4p without the density pressure 
function term (the DLSS model [8]), the positive classical solutions are proven locally in-time in 
one-dimensional periodic domain [2], and the global existence of a spatially periodic solution 
and its exponential convergence to an equilibrium state is shown [3] in terms of the entropy 
method and the Csiszar-Kullback inequality for "small" initial data. This is mainly due to the 
failure of the maximum principle which makes it impossible to establish a-priorily the upper 
and lower bounds of the density and obtain the global in-time existence of solutions with the 
strictly positive density. This, however, leads to the interesting results on the global existence 
of a nonnegative weak solution, which is first established in a one-dimensional bounded domain 
with the Dirichlet and Neumann boundary condition 14] where an interesting entropy estimate 
is introduced to show the global existence. Since then, some additional first order entropies 
are also obtained [161 [2T] . More recently, the global existence of multi-dimensional nonnegative 
weak solutions and their exponential decay to an equilibrium state is also shown for the DLSS 
model in a periodical domain [17J based on the extended multi-dimensional algorithmic entropy 
construction argument, and for the DLSS model with an additional term of the given drift 
potential [11] in the framework of variation and Wasserstein's metric subject to the finite initial 
mass. For more analysis related to the DLSS model or Eq. (II. 4p about numerical simulations or 
long time convergences, one can refer to the recent papers [U [T2| [TBI [23] and references therein. 
As for the bipolar quantum model (|1.2p . the existence of a stationary state is only analyzed 
recently [26]. Some interesting quasi-neutral limit has been analyzed recently [22j . 

However, there are few results on the global existence of classical (strong) solutions with 
the strictly positive density and the long time asymptotical behaviors of classical solutions for 
the macroscopic quantum models (II. 2|) and ()1.4|] in the whole spatial space, although there 
are a short time classical solution with the positive density for the DLSS model [2] and a 
global existence for nonnegative weak solutions [14^ IT6| |2T] ITU I18j . The main difficulties in 
dealing with the macroscopic quantum models (II. 2p and (II. 4p consist of the strong nonlinearity. 
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the degeneracy at vacuums, and the failure of the maximum principle, and the coupling and 
interaction between the two carriers for the bipolar case. 

It should be noted that it seems not obvious how to generalize the framework of entropy 
estimates and/or Wasserstein's metric, used for instance in [TTl [TU (TB] to establish the global 
existence of nonnegative solutions with the finite initial mass, to show the global existence of 
a classical solution with the strictly positive density for Eqs. (jl.2p or (jl.4p in the whole spatial 
space subject to the infinite initial mass (the case to be dealt with in the present case), since in 
general Poincare's inequality failed and it is not trivial to establish the a-priori uniform control 
of the density with respect to the time in order to understand the long time behavior of global 
solutions. 

In this paper, we are interested in the large time asymptotical behavior of the solution to 
the IVP ()1.2p - ()1.3p for the bipolar case and the IVP ()1.4p -( fL5]) for the unipolar case, and we 
shall show that the global classical solutions to IVP (fOjl - ffO]) and IVP (fLl]) - (fT3]) admit the 
character of self-similarity in large time so long as it is around the self-similar wave initially. To 
be more precise, let's introduce the quasi-linear parabolic equation 

Pt=p{p)xx, p'{p)>0. (1.6) 

It is well-known that Eq. ()1.6p has a unique self-similar solution W{x,t) up to a position shift 
(seel) 

pix,t)=:W{0, e = 7^, W{±^) = p±, (1.7) 

and the solution W{S,) is increasing if /9_ < and decreasing if p_ > p+, and — > as 
Id ^oo. 

To begin with, let us consider the case of regular initial data and take the initial datum of 
the IVP (fr2l) ~ (fL3|) and IVP (fLill -(fT3|) close to the self-similar solution W in some Sobolev 
norm. Then, we show in terms of the energy method below that there exists a global unique 
classical solution to the IVP (|1.2p - (|1.3p or the IVP (|1.4p ~ p.5p . which in particular tends to the 
self-similar wave W in large time with an algebraic decay rate. 

We should also mention that it is not obvious so far how to generalize the entropy functional 
or Wasserstein's metric technique as used in [HI [I2l [HI [23] to prove the convergence of global 
nonnegative solutions of the Eqs. (jl.2p and (jl.4p to the self-similar wave since the self-similar 
wave W is not a solution of the fourth-order quantum models (|1.2p and p.4p and the convergence 
itself is a singular process in large time. Moreover, some additional information on the lower and 
upper bounds of the density are needed to prevent the possible appearance of the singularity of 
the fourth order differential operator near vacuums. This is nontrivial however due to the failure 
of the maximum principal theory for the fourth-order equation. It is also interesting to whether 
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or not the global classical solution of the quantum models (jl.2p and ()1.4p shall converge to the 
self-similar wave W in large time for general initial data, instead of the small perturbation of 
the self-similar wave, it is left for further investigation. 

We first investigate the long time asymptotical behavior of global solutions to the IVP for 
the unipolar equation (ll.4p ~ ()1.5p . and then discuss the corresponding IVP for the bipolar model 



Let 



then 



where xq is determined by 



zo{x) = {po{y) -W{y + XQ))dy, 

J — oo 

zqx = Po{x) - W{x + Xq), 



+00 

{po{x) - W{x + xo))dx = 0. 

-00 



For the unipolar case, our main result on the global solution and its large time behavior of 
IVP (fLl]) - (fT3]) is given as follows. 

Theorem 1.1 Let p'{p) > for p > 0. Assume that 6 =: |p+ — p-| ^ 1 and zq e H'^{R) 
with 5o =: \\zQ{x)\\jj:i^-^-j small enough, then there is a unique global strong solution p > of 
IVP ([LlD-dLS]) such that 

p-W e L°°([0, 00); H'^{R)) n l2([0, oo); H^{R)), 

and the solution p converges to the self-similar wave W{^^=) of the Eq. (jl.6p with an algebraic 
time decay rate 

\\d^,ip - wmwiHR) < c{i + 1)-"^, k = o,i, 2, 

||(p-PF)(t)||ioo(K) ^C(l + t)-i, 
where C is a positive constant dependent of 5 and 60 . 

Next, we state the main result on the convergence to the self-similar wave for the bipolar 
case. Although it leads to additional difficulties, the coupling between carriers in l\l.'2\] may 
cause some cancelation, and it is not clear that both the densities of the bipolar QDD (II. 2|) 
behave still or not like the unipolar one for a small perturbation. 

Denotes 

zi{x)= I {pifl{y)-W{y + XQ))dy, i = a,b, (1.9) 
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then 



where xq is determined by 



Zqx = Pi,oi^) -^(^ + ^o)^ i = a,b, 



(Pifiix) - W{x + xo))dx = 0. 



Theorem 1.2 Let p'{p) > for p > 0, and 5 =: |p+ — P-\ ^ 1- Assume that mixeRPi,o > 

{i = a,b) with f{pa,o — Pb,o)dx = 0, and Zq G ff^(M) with 6o =■ \\zq\\h3(r) + ||-2o||//3(r) small 

enough, then there is a unique global strong solution {pa,Pb,E) of the IVP (|1.2|) - (|1.3|) with 
Pa > 0, Pb > such that 

Pi-W e L°°([0, oo); H^{R)) n L^{[0, oo); //^(M)), i = a,b, 

and both pa and p^ converge to the self-similar wave of the Eq. (II. 6p with an algebraic 

time decay rate 

\\dl{pa - W)(t)||L2(K) + \\dl{p, - H^)(t)||i2(M) ^ C(l + t)-^, = 0, 1, 2, 

\\{Pa - W^)(t)||L-(M) + \\{Pb - W^)(t)||L-(R) < C(l + t)-i, 

mmmim'^ce-^', 



where C > and /? > are constants dependent of 5 and 5o . 

Notations. L^(M) and H^{W) denote the usual Lebesgue integrable functions space and 
the Sobolev space with norm || • ||lp(]k) and || • ||/^fc(K) respectively, we also use || • || to denote 
II • ||l2(]8-) for simplicity. C and c are used to denote general positive constants. 

2 Proof of main results 

We shall prove Theorems 1 1 . HfL2] in this section. The key is to establish the a-priori estimates for 
short time strong solutions. Without the loss of generality, we establish the expected estimates 
in order to prove Theorem 11.11 in Sect. 12. H and show how to derive the estimates about the 
electric filed in the proof of Theorem 11.21 in Sect. 12.21 



2.1 The unipolar case 

In this section, we shall transform the primary equations in order to study the existence and in 
particular its large time behavior of the global solutions of the IVP ()1.4|) - (ll.5p . 
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Denote 



then 



.(.,0 = /_Jp(.,«)-H--(^))*. (2^1) 

z, = p{x,t)-W{!^^). (2.2) 
We will derive the fourth order parabolic equation for z. Since W satisfies 

Wt = p{WU, (2.3) 
and p = p{x,t) satisfies Eq. ()1.4p . we have 

{p - W)t - {p{p) - p{W)),, + (p f ) = 0. (2.4) 



Integrating ()2.4p over (— oo, x) with respect to the spatial variable and assuming pxx ^ as 
|x| oo, we obtain from (j2.ip and ()2.2p the parabolic equation of the fourth order for z of the 
following form 



with the initial datum 



where 



£2 

Zt - {p{W)Zx)x + -Zxxxx = ifl + f2)x, (2.5) 



z{x,0) = zo{x), (2.6) 



h=p{^x^W)-pl[W)zx-p{W). (2. 
Note that we have used the fact 



The existence of the global solution and the large time behavior for the IVP (|2.5p - (j2.6p is 
obtained by the following proposition. 



Proposition 2.1 Let p'{p) > for p > and 6 = |p+ — P- \ <C 1. Assume that \\zo\\h3(k.) ^ So 
with 6o > sufficiently small, then there is a unique global strong solution z to the IVP ()2.5p - 
(j2.6p satisfying 



z G L°°([0,oo);/7^(M))nL2([0,oo);F^( 
||9^z(-,t)||i2(K) ^C7(5,<5o)(l + t)"i A: = 0,1,2,3, 
\\d^,zi;t)\\Lo^iR)^C{6,6o){l + t)-"-^, k = 0,1,2, 

where C{5,6o) > is a positive constant depending only on 6 and 5q. 
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Remark 2.2 It is sufficient to prove Theorem \l.l\ in terms of Proposition \KT\ due to the relation 
between p and z 

p = W + 

and the transformation of i2.S^) - [275\) . The positivity of p can be assured by the positivity of W 
and the smallness of z^- From the above, we also have 

||a^(p-Ty)(t)||i2(K) = ||a,^+iz(t)||i2(K), fc = o,i,2. 

In order to prove Proposition [2?T1 let us assume that for the local in-time solution and T > 

5t= max sup (1 + t)t ||a^z|| < 1. (2.9) 

fc=0,l,2,3o^t^T 

By the Nirenberg's inequality and the above assumption, we have 

l|9^IUoo(K) <c5T(l + t)-^, fc = 0,l,2. (2.10) 
The theorem for the existence of the local in-time solution is 

Theorem 2.3 Letp'(p) > for p> and assume that |Uo||_ff3(']R\ small enough and inf (zq^ -|- 
W{x + xq)) > {i.e. inf po > 0). Then there exists a T* > such that there is a unique local 
solution of the IVP ([23]) - (p:6]) for t £ {0,T*) satisfying 

\\zi-,t)\\m{R) < oo, p{x,t) = z^{x,t) + W{x + Xo,t) > 0. 

The proof of Theorem 12.31 can be obtained by a standard method (see |llj). we omit it. 
Now, we list the L^-estimates of the derivatives of W and n = \/W as follows. 



Lemma 2.4 Let W be the self-similar solution of il.b]) and n = \/W , then it holds that (seefBj) 

II W-,i)llLP(M) ^ C5(l +t)~2+i, (2.11) 
||5X-,t)||z,P(R) ^C<5(l + t)-^+^, (2.12) 

for j > 0,p E [l,-|-oo], C > is some constant. 

For /i, /2, we have the following estimates. 

Lemma 2.5 Under the assumption ^2. 9\) . it holds for /i,/2 



fi = 0{6T + S)z,, + 0{6)r2, (2.13) 
/2 = 0{6t)z^, (2.14) 
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where the function rj.[x,t) is related to the kth order derivative of W with respect to x, which 
satisfies by the definition 

lkfc(-,t)||LP(M) ^ fc = 0,l,2,---, andpG [l,+oo]. (2.15) 

Proof. From Lemma |2.4| (j2.7p and (j2.8p . we have 

•'^ 2 W + z^ 2 """^ 

= 0{l)iW^ + 2W,-z,, + zl, + W,,) 

= {2W^ + z^^) ■ z^^ + 0{5)r2 

= 0{5t + 6) ■ z^^ + 0{5)r2, 

f2 = p{z, + W)-p'{W)z,-p{W) 

zl + 0{zl) 



2 

= Cz^{zj, + z^) = 0[bT)z^, 
by which the proof is easy. □ 
Lemma 2.6 Under the assumption \2.y\) . it holds for the local in-time solution z 

\W)f + r WzMfds + r ^ 0(5 + 5o)', (2.16) 



for ^ t ^ T , provided that 5t + 5 is small enough. 



Proof. Taking the L^-inner product of (j2.5|) with z, we get with the help of the integration 
by parts 

/ [zt - {j>iy^)z^)x + —ZxxxxA ■ zdx = Zf z + p'{W)zl + -T-zl^dx 

= ^^Ikf + ylkxxf + ^p'(W^)^'dx. 
By Lemma [231 for /i,/2 and Cauchy's inequahty, we get 
/ {fi + f2)x-zdx = / -{fi + f2) ■ Zxdx 

^ a\\zxf + Om\hf + \\hf) 

^ a\\zxf + 0{5t + 5f{\\zxxf + \\zxf) + 0{5^){l + t)^i, 

where a > is a constant such that a + 0{5t + Si)"^ ^ tti inf p'iW). Combining these estimates, 
we get 

\^M? + - + 5f)\\zxx\? + I {p{W) -a- 0{6t + 5f)zldx ^ 0{6''){l + *)"§, 

I at I Its 
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Since 6 + 6t <. 1, we have — - 0{S + 5t) = O(e^) and 

^ - zf + 0{e'')\\z^^f + [ ip'iW)-a-0{6T + 5))zldxi:0{5^){l + t)-l. 



Integrating the above inequahty with respect to the time from to t, we get 



\z[t]\?+ / \\z^^{sy^' 





''ds+ f I p{W)zl dxds ^ O(<5o + 5f. 
Jo Jr 

This, together with the fact p'{W) > for > 0, gives (f2T6]) . □ 
Lemma 2.7 Under the assumption ^2. 9\) . it holds for the local in-time solution z 

+ r Wz^Mfds + r Wz^^Mfds ^ 0(5 + <5o)^ (2.17) 

JO JO 

(1 + t)\\z,{t)f + / (1 + s)(||z,,(s)f + \\z^Us)f)ds ^ 0(5 + 5o)2 (2.18) 

JO 

/or ^ t ^ T provided that 5t + § is small enough. 

Proof. Differentiating the equation (12. 5p with respect to x and taking the L^-inner product 
of the resulting equation with z^ , we get in view of the integration by parts 

(/i + f2)xxZxdx 

[{zt)x-{{p{W)z^)^)^ + —{z 

xxxx)x\ * ^x^"^ 



-- — \\z + — lb IP 



+ f p'{W)zljx + h, (2.19) 



where 



Kil = I / p" {W)W^z^Zr.dx\ ^ 0{l)\\W^\\loo\\z,f + a\\z,^f , (2.20) 

JM 

with the same a as in Lemma 12.61 By Lemma 12.61 and Cauchy's inequahty, we get 

I / (/l + f2)xxZxdx\ =1 / (/i + /2)a.Za;a;(ix| 
JM JR 

<.0{l)\\{fi+f2)xf + a\\Zxxf 

^0{5t + 5)(||2xxf + 11^0=0=0=11') + a\\z^xf + 0(52)(1 + t)-f , (2.21) 

where we have used the assumption p.lOj) and a similar analysis as in Lemma [2. 51 with the help 
of Lemma 12.41 Combining (|2.19p - (|2.21|) and integrating the resulting inequality over [0, t\ , we 
have after simplifying 



t rt 

2' / '\z^^{s)fds+ \\z^^^{s)fds^O{5o + 5f, 

JO 
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which gives (j2.17p . To obtain (j2.18p . differentiating the equation (j2.5p with respect to x and 
taking the L^-inner product of the resulting equation with (1 + s)zx, similarly with the former 
we get after the integration by parts 

/ (/l + f2)xxZx{'^ + t)dx = / [{Zt)x - {{p {W)Zx)x),j, + %-{Zxxxx)x\ ' Zx{l + t)dx 
JR JR ^ 

=\^y{X^l)\\zxf\--\\zxf 

+ -(1 + t)\\zxxxf + / (1 + t)p'{W)zlxdx + (1 + t)h, (2.22) 
^ Jm. 

with the estimate 

|(l + t)/i|=| / {l + t)p"{W)WxZxZxxdx\ 

^0(1)(1 + mWxWlo. \\zxf + a(l + t)\\zxxf 

^0{l){6^)\\zxf + a(l + t)\\zxxf, (2.23) 

where we have used Lemma for ||Wa;|||^oo ^ C(5^(l + t)^^. We also have estimates 

/ (1 + i)(/i + f2)xxZxdx =1 / (1 + t)(/i + f2)xZxxdx\ 
JR JR 

^0{5t + -5)(1 + Odl^xxf + ll^xxxf ) + a(l + OIkxxf 
+ 0(52)(l + t)-i. (2.24) 

Combining (j2.22p - (|2.24p and integrating the resulting inequality over [0, ti] will give us after 
simplification 

rh 

(1 + h)\\zx{ti)f + / (1 + m\zxx{t)f + \\Zxxx{t)f)dt ^ 0{5 + 6of, 
JO 

where we have applied the derived estimates in Lemma 12.61 to get that 

r \\zxfdt^O{6 + 6of 

JO 

for ti G [0, T]. This gives (I2J8D . □ 
Based on Lemma 12.61 and Lemma \'2.7\ we can perform the estimates of the second and the 
third order derivatives of z in the same procedure and we have 

Lemma 2.8 Under the assumption ()2.9p . it holds for the local in-time solution z 

(1 + tf\\Zxx{t)f + [\l + sf{\\zxxx{s)f + \\Zxxxx{s)f)ds ^ 0{S + <5o)', 
JO 

(1 + tf\\Zxxxit)f + / (1 + sfi\\Zxxxxis)f + \\Zxxxxxis)f)ds ^ 0{6 + 5of 

Jo 

for ^ t ^ T provided that 5^ + S is small enough. 
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Proof. We give the sketch of the proof. Performing J^^ f^{l + 1)"^ (2.5) xxZxxdxdt, we can get 
as in Lemma 12.71 

rh 

{l+tlf\\Zxx{t)f + / {l + tf{\\Zxxx{s)\\^ + \\Zxxxx{s)\\^)dt 

Jo 

^ / {l+tf\\zx\\'i^\\zxx\?'dt+ / {l+tf\\rifdt, (2.25) 



with r4 defined in Lemma 12.51 satisfying 

By ()2.10p we also have H^^xHioo ^ C(5|,(l + 1)^"^. Thus, by (|2.25p and Lemma [2?71 we have 

{l + tif\\zxx{t)f + / {l+tf{\\zxxx{s)f + \\zxxxx{s)f)dti^O{5 + 6of. (2.26) 

JO 

Similarly, by performing J^^ /r(1 + i)^ {'^■^)xxxZxxxdxdt and with the help of all the derived 
a-priori estimates we can obtain 

{l+tf\\Zxxx{t)f + [ {I + sf{\\Zxxxx{s)f + \\Zxxxxx{s)f)ds ^ 0{6 + 6of. 

Jo 

□ 

The proof of Proposition \2.1\ The Lemmas 12 . 6H2 . 8 1 show that the local solution satisfies the 
uniform bounds for short time (5 <C l,5o ^ 1) when the initial perturbation is small enough. 
By using continuous methods, we can extend the local solution to be a global one, which also 
satisfies Lemmas I2.6H2.8I for any time. Then the existence of the global solution to the IVP 
(j2.5|) - (|2.6p is proven. The Theorem II. II is a direct corollary of Proposition 12.11 □ 



2.2 The bipolar case 

In this subsection, we prove Theorem 11.21 about the IVP (ll.2p - (ll.3p for the bipolar case. We 
also shall transform the primary equations in order to study the existence and in particular its 
large time behavior of global solutions of the IVP (|1.2p - p.3p . 
Denote 

z\x,t) = £ (^p^{y,t)-w(^^^^^^^dy, i = a,b, (2.27) 



which implies 



zl = p,ix,t)-W{^^], i = a,h. (2.28) 
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Similar to the unipolar case, assuming {pi)xx — > as |x| — > oo, we can rewrite the equations 
(fr2D as 



^ zi - {p'{W)zl)x + y4x.. + (-l)*+'(4 + W)E = (/o + /,,2)., ^ = a,b, 

Et - {p'{W)Ex)x + ^-^Exxxx + 2WE = {fa,l - h,l + fa,2 - /b,2)x " (4 + 4)E, 



(2.29) 



with initial data 



where 



z\x, 0) = zlix), 4(±oo) = 0' i = a,b, (2.30) 



{pa{y,t) - pb{y,t))dy = z" - z\ 

-oo 

' 2 W + zl, 2 ' 
/i,2 =P(4 + ^) - P'{W)zi - p{W), i = a,b, 

and we recall we make use of the symbols (—1)° = —1 and (—1)'' = 1 for the simplicity of 
statements. 

We have the following results: 

Proposition 2.9 Let p'{p) > for p > and 6 = |p+ — P-\ ^1- Assume that ||4(^)ll-f/3(R) + 
ll4(-^)llii'^(lR) ^ ^0 with 6o > sufficiently small, then there is a unique global strong solution 
{z'',z'',E) to the IVP dOoD - dOn]) satisfying 

(z^z^) G (L-([0,cx));//3(M))nL2([0,cx));/75(M)))2, 
||9^(z^ z')(-, t)||i2(K) ^ C{6, 6o){l + t)--2, k = 0,l, 2, 3, 
||a^(z^ /)(., t)llioo(M) ^ C{5, 5o){l + t)-"^, k = 0,l, 2, 
\\E{t)U2^u)^C{5,6o)e-^\ 

where C{5,6o) > and (3 are some positive constants depending only on 5 and 5q. 

The proof of Proposition 12.91 can be made in the similar fashion as Proposition 12.11 Here, 
we only show how to deal with the electric field and show its exponential decay. Let us assume 
that it holds for local in-time solutions that 



5t= max sup (1 + i)^ < 1, for T > 0. (2.31) 

fe=0,l,2,3 o^4^T 
I— a,D ^ ^ 



Similar to Lemma 12.51 we have the following properties for fi^i and /j^2) whose proof is very 
similar to that of Lemma 12.51 we omit details. 
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Lemma 2.10 Under the assumption (j2.3ip . we have for i = a,b and j = 0, 1, 2, that 

difi,i = 0{5t + 8)di-'^z' + 0{5)r,+2, difi,2 = 0{5T)di-'^z\ 

di{fa,l - fb,l) = 0{5t + d){dt'E + dt^E), di{fa,2 - h,2) = 0{dT)di+'E, 



where the function rf^{x,t) is the same one as in LemmaW. 



We also have the bipolar version of Lemmas I2.4H2.8I as similar ways as in the unipolar case. 
Indeed, we can obtain the following a-priori estimates. 

Lemma 2.11 Under the assumption (j2.3ip . it holds for the strong solutions {z"',z^,E) that 



dl{z\z'){t) \jzf dl^\z\z'){s, 



ds + 



* 2 

dlE{s) dsi^O{6 + 5of, 







(2.32) 
(2.33) 



for k = 0,1, 2, 3, and 



d^E{t) ^ 0{6 + (5o)e"^'=*, fc = 0, 1, 2, for some constants Po> (3i,p2> 0, 
for t G [0, T], provided that 6t + 5 is small enough. 



(2.34) 



Proof. We only give a sketch of the proof here since it is very similar to the unipolar case. 
Taking the kth order derivative of the first equation in (j2.29p with respect to the spatial variable 
X, and taking the inner product of the resulting equation with d^z^, we get for k = 0,1, 2, 3, 
with the help of ([23T]) and Lemma [13] and [230l that 

.2 



E 



1 d 
2di 



+ E / p'iwM+'z^dx + '-Y^ 

I,— u,,u i=a,b i=a,b 

(_l)5o. V / ^'i~^^^^-'^^^\f^^ + f^,2)^t^-'"'^'dx 
t=a,bJ^ 

k-1 .. 

-{l-5ok) E E^U dt-'-\p"{W)WM-'^z^dl+^z'dx 



+ / Wid^Efdx 



i=a,b j=Q 



1 



+ ^i-l)HSok + Sik) I {z" - z')Ed',+'+'z-dx 



{S2k + Ssk) I dt~\{zl-zl)E)d':,+^z^dx 

k-l 



f di^^{zlE)df-^--^''Edx - S^Ci [ d'^-^WdiEdlEdx 
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+ 0{6T + Sf H 



,2+max(fc-l,0) i 



i=a,b 



i=a,b 



k 



k-l 



max{k, 1) 



i=a,b j=0 

+ {S2k + Ssk) [0{d^T){\\Ef + 63k \\E^f) + as 

r, k-l 

+ 04 



j=0 



where 6ij = 1 if z = j and Jjj = if i 7^ j, Oj's are small positive constants such that for 
sufficiently small 6t, S and T , p'{w) — a\ — — 0{6t + (^)^ > 0, ^ — as — 0{6'^) > and 
W-a4- 0(4) > 0. Thus, it is easy to obtain the desired result ()2.32p by considering, in turn, 
every k = 0,1, 2, 3, respectively. 

Now, we prove the second result (12.340 . Taking the fcth order derivative of the second 
equation in ()2.29p with respect to the spatial variable x, and taking the inner product of the 
resulting equation with d^E, we have for A; = 0, 1, 2, in view of (I2.3ip and Lemmas 12.41 and 12. 101 
that 



1 d_ 

2di 



dlE{t) + Jofc / p{W)Eldx + - 



+ 2 / Wid'iEfdx 



=6ik I p'{W)E,E,^,dx + 62k [ p"{W)W^E,E,,,,,dx - 2(1 - 6ok){-l)'' [ WEdfEdx 
Jr Jr Jr 

+ 25ik [ W{d^Efdx - 252k [ WE^d^Edx + 25ok [ (z" + z^)EE^dx 
Jr Jr Jr 

-Sik / {z" + z^)(E^E^^ + EE^^^)dx + 52k / (z" + z^){2,E^^E^^^ + E^E^^^^)dx 
Jr Jr 

- (-1)'=(1 - 52k) [ (/a,l - fb,l + fa,2 " fb,2)df+^Edx 

Jr 

+ hk j {fa,l - fb,l + fa,2 - fb,2)xd^Edx 

Jr 

^5ife[0(l) \\E^f + as WE^xxxf] + 52k[0{6^) \\E,f + ag \\E,xxxf] 



+ {l-6ok)[Oi5^)\\Ef + a-r 



dV'E 



+ 52kO{5^)[\\E£ + \\E^^^\ 



+ 5ok[0{5l) WeW" + as WE^f] + 5ifcO(4)[||S,f + \\E^^f] + 5ik[0{5l) + og \\E^ 
+ hkO{5^T)[\\Exx\? + \\E^xx\?] + 52fc[0(4) \\Exf + aio \\E^xxx\?] 
+ 5ok[0{5 + 5t) ll^.f + 0{5 + 5Tf ll^.f + an \\E,xf] 
+ 5ik[0{5 + 5Tf{\\E,\\' + \\eU\) + ai2 \\E,xx\?] 
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+ 52k[Oi5 + 5T)H\\E..f + + ai2 \\E,,,,\\% 

for some sufficiently small constants aj's. Thus, by considering every = 0, 1,2 in turn, there 
exist /3o > /3i , /32 > for 6 and T small enough such that 



1 d 
2di 



^{l-5Qk)0{5T + 5)\\Ef , fc = 0,l,2. 



From Gronwall inequality, we can obtain 

d^^E{t) < 0((5o + (5)e-^^*, A; = 0,1, 2. 

Finally, let us prove the last result ()2.33p . Taking the A;th order derivative of the first equation 
in (j2.29p with respect to the spatial variable x, and taking the LP' inner product of the resulting 
equation with (1 + t)^d^z\ we get for k = 1,2,3, with the help of i^Mil, (1Z32]), (fOil) and 
Lemma El and Eini that 



d 
di 



(1 + 0'' 



i=a,b 
2 



,k+2 i\ 
X ^ 



i=a,b 
2 



(1 + *)' Y.Y.Cil dl-h'{W)di-^'z^dt^'z'dx 

.-u,,„ i=a,b j=0 ''^ 

2{l + t)' [ W{d',Efdx-{l+t)'^ f dt'izl + W)diEd',Edx 

■^^ i=a,b j=0 •^'^ 

{l + tfY. i zl{d';Efdx + {l + tfY. I dtHfl + f2)d'x'-'z^dx 



k 



+ E \H 



2 

i=a,b i=a,b 

+ S2k{i + tf [o^-^') II4.II' + «i5 + o{i) I 

i=a,b 

^xx 1 1 1 1 ^xxx 1 1 I ( ) II ^xxxx 1 1 ^16 1 1 ^xxxxx 

i=a,b 

+ 0{6 + 6T)Oi6 + 6of, 
which implies the desired result as the way as the previous. 



□ 



The proof of Theorem and Proposition \2.9l Since the proof is very similar to the unipo- 
lar case in view of Lemma l2.11l we omit details. □ 
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